Determining the in-plane orientation of the ground-state orbital of CeCu 2 Si2 
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We have successfully determined the hitherto unknown sign of the B\ Stevens crystal-field param- 
eter of the tetragonal heavy- fermion compound CeCu2Si2 using vector q dependent non-resonant 
inelastic x-ray scattering (NIXS) experiments at the cerium N4,5 edge. The observed difference 
between the two different directions q||[100] and q||[110] is due to the anisotropy of the crystal-field 
ground state in the (001) plane and is observable only because of the utilization of higher than 
dipole transitions possible in NIXS. This approach allows us to go beyond the specific limitations of 
dc magnetic susceptibility, inelastic neutron scattering, and soft x-ray spectroscopy, and provides us 
with a reliable information about the orbital state of the 4/ electrons relevant for the quantitative 
modeling of the quasi-particles and their interactions in heavy-fermion systems. 

PACS numbers: 71.27.+a, 75.10.Dg, 78.70.Ck 



4/ heavy-fermion systems exhibit many fascinating 
phenomena including unconventional superconductivity 
and non- Fermi liquid behavior emerging in the vicinity of 
quantum critical points. In these systems the conduction 
band electrons are dressed by the hybridization with 4/ 
electrons to form Landau quasi-particles with large effec- 
tive masses. It is of great interest to identify the relevant 
quantum state of these 4/ electrons in order to under- 
stand the formation and character of these quasi-particles 
and their interactions. Specifically, it is well known that 
the crystal-field split Hund's rule ground state of the 
unhybridized 4/ electrons is highly anisotropic and it 
was recognized already early on that this anisotropy is 
important for the modeling of the properties of heavy- 
fermion systems [l], H| ■ There are examples for strongly 
anisotropic heavy fermi liquid formation Q and several 
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FIG. 1: (color online) Angular distribution of two Ft crystal 
field states that have the same | ±5/2) admixture (a 2 = 0.81) 
but different signs of a. 



authors have pointed out how the crystal-field anisotropy 
affects the band structure and shape of the Fermi sur- 
face |4[ and / or how it mayinfiuence the formation of the 
superconducting state [5j-l7|. 

However, for the largest group of heavy-fermion ma- 
terials where the Ce 3+ (or Yb 3+ ) ion has tetragonal 
point symmetry (e.g. all CeM2Si2, CeM2Ge2, CeMIris, 
YbM2Si2, YbM2Ge2 with M being a transition metal, 
or CePtaSi) the crystal- field ground-state is only incom- 
pletely determined because the orientation of the ground- 
state orbital in a symmetry higher than twofold cannot be 
determined with standard dipole techniques like inelastic 
neutron scattering, or more recently applied for crystal- 
field purposes, soft x-ray absorption spectroscopy. In this 
work we therefore utilize the higher than dipole transi- 
tions in non-resonant inelastic x-ray scattering (NIXS) 
and introduce it as a spectroscopic tool new to this 
research field of heavy fermion physics. We will show be- 
low that we can resolve reliably the so-far unknown orien- 
tation of the ground-state orbital in the archetype system 
CeCu2Si2 which was the first heavy-fermion compound 
where unconventional superconductivity was discovered 
more than 30 years ago (l9(. 

In tetragonal 4/ systems the sixfold degenerate Hund's 
rule ground state of Ce 3+ (J=5/2) is split into three 
Kramer's doublets under the influence of the crystal 
field. The crystal-field Hamiltonian can be written in 
Stevens approximation as H C f = B%0% +B\0% + B\0\ 
and the eigenfunctions can be represented in the ba- 
sis of \ J Z ) when the fourfold symmetric tetragonal [001] 
axis is chosen as quantization axis. There are two Tj 
d oublets T\ = a\± 5/2) + y/i - a 2 \ qp 3/2) and T 2 - = 
y/l — a 2 1 ± 5/2) — a\ =F 3/2), and one Tq which is a 
pure | ± 1/2) doublet. The as a pure | ± 1/2) has 
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full rotational symmetry around [001] but the mixed T7 
states do not. Both T7 states have a fourfold symme- 
try around [001] and for a given spatial distribution of 
the 4/ wave function there are two solutions which differ 
in their orientations within the (001) plane by 45° (see 
Fig. 1). Which orientation applies to the ground state 
depends on the sign of a. For a > the wings of a T? 
ground state point along [100] and for a < along [110]. 
The three Stevens parameters B®, B®, and B\ are fully 
determined by the mixing parameter a and the energy 
level splittings AEi and AE2, whereby the sign of a is 
related to the sign of B\ such that a < corresponds to 
B\ > and vice versa. 

For CeCu2Si2 the crystal-field transition energies have 
been determined with inelastic neutron scattering [2(3] 
on polycrystalline samples. Inelastic neutron scattering 
is the ideal tool to determine the level splittings but 
the combination of phonon scattering and broad mag- 
netic excitations in the same energy window prevents 
the unambiguous determination of the magnetic inten- 
sities. Therefore the ground-state wave function was de- 
termined with single crystal susceptibility measurements 
20l . 2~fl ] and the crystal-field parameters which reproduce 
the susceptibility are given by Goremychkin et al. [13] 
and correspond to |a| = 0.88. Note, that the sign of a 
was not determined. Single crystal susceptibility mea- 
surements are a common way of determining wave func- 
tions in heavy fermion compounds. More recently, soft 
x-ray absorption spectroscopy with linear polarized light 
was shown to be a useful local probe to determine the 
anisotropy of the wave functions spectroscopically 22- 
Eri ]. But both soft x-ray absorption spectroscopy and 
neutron scattering are governed by dipole transitions so 
that these methods are insensitive to anisotropies with 
a higher than twofold rotational axis and therefore can- 
not distinguish the two possible orientations of the T-? 
orbital in the tetragonal (001) plane, which is synonym 
with not being able to determine the sign of a or B\. 
In principle this in-plane anisotropy can be observed in 
direction dependent isothermal magnetization measure- 
ments, but for most cerium materials it would require 
even at 1 K magnetic fields of more than 20 T. Therefore, 
as for many other tetragonal heavy- fermions compounds, 
the ground state in the archetype heavy-fermion com- 
pound CeCu2Si2 could only be characterized by the abso- 
lute value of the mixing parameter |a| (or \Bf\) leaving 
the question about the orientation of the ground state 
within the (001) plane unanswered. 

Here we present NIXS data of CeCu2Si2 single crystals 
and show how the higher multipole transitions can reveal 
the anisotropies of the ground-state wave-function which 
have not been accessible with other spectroscopies. NIXS 
has became feasible thanks to the high brilliance of mod- 
ern synchrotrons and advanced instrumentation and has 
developed to a rapidly growing field providing comple- 
mentary means to x-ray absorption methods in the study 
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FIG. 2: (color online) Left: k th order term of the radial 
part of S(q,uj) as function of momentum transfer. Right: 
k*' 1 order contribution of angular part of the scattering func- 
tion expressed in terms of S(cj) versus energy transfer for 
(R f \j k (qr)\Ri) = 1 (see Eq.rj). 



of core level excitations, but offering the possibility to go 
beyond the dipole limits [§-[3 • For rare earth materials 
higher multipole Ad — » 4/ core level excitations (N^- 
edge) have been observed for large momentum transfers 
|q| by Gordon et al. [Io1.[lr3] and Bradley et al. [Hj]. It was 
actually demonstrated experimentally how the dipole sig- 
nal decreases with increasing |q| while dipole forbidden 
transitions appear. These findings are in good agreement 
with simulations based on the local many body approach 
by Haverkort et al. [9( and it was suggested that vector 
q dependent NIXS experiments with large momentum 
transfers on single crystals should have the potential to 
give insight into the ground state symmetry. The direc- 
tion of the q vector with respect to a crystallographic 
axis is the analog of the polarization in an x-ray absorp- 
tion experiment with linear polarized light. Gordon et 
al. performed first NIXS experiments on cubic single 
crystals of MnO and Ce0 2 at the Mn M 2 , 3 and Ce N 4;5 
edges and found some differences for different q direc- 
tions in the higher multipole scattering. This has been 
interpreted as sensitivity to anisotropies in charge densi- 
ties with higher than two- fold symmetry [111 ]. 

For the sake of clarity, we first review briefly the theo- 
retical des crip tion of inelastic x-ray scattering 0, [l(J 13, 
llEl [Lsl I27M29I] and show later how the crystal-field affects 
the scattering cross section. The double differential cross 
section is the product of the Thomson photon cross sec- 
tion (^§i) Tho and the dynamical structure factor S(q,ui), 
the latter containing the physics of the material under 
investigation: 
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The dynamical structure factor is a function of the scat- 
tering vector q = k; — kf and the energy loss uj = u>i — uif 

S(q,w) = ]T I </|e 4qr |*) -tof- M- 
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FIG. 3: (color online) Simulations: Comparison of in-plane 
and out-of-plane scattering function S(q, oj) of pure J z states, 
in-plane q||[100] (blue), out-of-plane q||[001] (red), and in 
between q||[101] (orange). The calculations are done for 
Iql =9.3 A -1 and are convoluted with a Lorenzian with 
FWHM = 0.3 eV and a Gaussian with FWHM = 1.32 eV. 
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FIG. 4: (color online) Simulation (top): Comparison of the 
scattering function S(q, uj) for two in-plane directions q||[100] 
(blue) and q||[110] (green) assuming a F7 ground state with 
a < 0. The calculations are convoluted with a Lorenzian with 
FWHM = 0.3 eV and a Gaussian with FWHM = 1.32 eV. Ex- 
periment (bottom): NIXS data of two CeCu2Si2 single crys- 
tals. The blue dots and lines corresponds to the set-up for 
q||[100] and the green ones to q|| [110] . The error bars reflect 
the statistical error. The lines are guides to the eye. 



Here i and / are the initial and final states. In order 
to distinguish between dipole, octupole etc. terms the 
transition operator e lqr is expanded in semi-normalized 
(Racah's normalization) spherical harmonics C^ m and 
C\ m . This results in a sum over spherical Bessel func- 
tions jk (qr) and the wave functions can be factorized into 
a radial and angular part so that S(q,w) can be written 
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Due to the triangle condition and parity selection rules 
only terms with k=l (dipole), 3 (octupole), and 5(triakon- 
tadipole) contribute to the Ce Ad — » 4/ (N^s) transitions. 
The radial part (jR/|jfc(qr)| Ri) of the wave functions have 
been calculated within the Hartree-Fock approximation 
using Cowan's code [13] and the k th contributions are 
shown as function of momentum transfer |q| in the left 
panel of Fig. [5] For moderate magnitudes of |q| the ra- 
dial part is dominated by dipole scattering, but already at 
5 A -1 octopole scattering is non- negligible, and at even 
higher momentum transfers the scattering is dominated 
by the higher multipoles. This behavior is commonly 
called q-dependent multipole selection rules. The right 
panel of Fig. [2] shows the k" 1 order of the angular part as 
function of energy loss. Higher multipoles have different 
selection rules so that extra intensity at different energy 
losses becomes visible in the angular part when at large 
|q| octupole and triakontadipole transitions take place. 

Having described the |q| dependence of the NIXS in- 
tensities, we now turn to the vector q dependence which 
is at the heart of our study. We first of all evaluate 
the sensitivity of NIXS at the N4.5 edge to crystal-field 
anisotropies in general by comparing simulated spec- 
tra for different directions of q. These crystal-field 
anisotropies are included when the initial and final states 
in Eq. [2] are eigenstates of a Hamiltonian that contains 
in addition to the atomic Coulomb and spin-orbit inter- 
actions also non-vanishing crystal-field terms. We note 
that the interference terms which vanish if the ang ular 
intensity is integrated over all directions q flil Il8j are 
included in our calculations. Figure [3] shows the simula- 
tion of S(q, lj) for large momentum transfers for the three 
pure J z states, considering some realistic lifetime broad- 
ening and width due to instrumental resolution (see be- 
low). The pure states have rotational symmetry so that 
we compare the in-plane scattering (q|| [100] ) with scat- 
tering out-of-plane, i.e. for q||[001], and some direction 
in between (q||[101]). There is a clear direction depen- 
dence, so that the different J z states are distinguishable. 
This is similar to the in-plane/out-of-plane polarization 
dependence in soft x-ray absorption at the cerium M4 5 
edge [22|,[23|. 

The next step is to investigate the in-plane anisotropy 
of a mixed Tj state as shown in Fig. [1] For that we 
calculate and compare S(q, uS) with a < for the two in- 
plane directions q|| [100] and q|| [110] and we end up with 
distinguishable spectra (see top of Fig. 0]). This clearly 
shows that such an experiment can be used to determine 
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the sign of a. The simulation in Fig.[4]assumes a negative 
value for a and the corresponding orbital orientation is 
shown in the inset. For a 45° rotation around the o 
axis [001], i.e. for a positive value of a the spectra are 
inverted. We found that these in-plane spectra are fairly 
insensitive to the precise value of a as long as a is neither 
zero nor one. In the latter case the orbital would have 
full rotational symmetry around [001] and consequently 
5(q, lu) would look the same for both in-plane directions. 

To determine the orbital orientation in CeCu2Si2 we 
performed Ce N^s NIXS measurements at the IXS end 
station of the Taiwan beamline BL12XU at SPring8 on 
two CeCu2Si2 single crystals. One crystal had a polished 
(001) and the other one a polihed (110) surface so that 
q||[100] and q||[110] could be realized in a specular ge- 
ometry. The crystals were cooled down to about 20 K in 
a closed cycle cooler so that only the ground state was 
populated. The incident photon energy of 9890 eV and a 
scattering angle of 28=138° correspond to a momentum 
transfer of |q| = 9.3 A" 1 . The energy resolution as de- 
termined from the full width at half maximum (FWHM) 
of the elastic line was 1.32 eV. The signals of nine spher- 
ically bent backscattering Si(555) analyzers at two meter 
distance from the sample were combined in order to gain 
higher count rates. The energy loss was determined from 
the energy position of the elastic line measured before 
and after each set of spectra. 

The measured Ce spectra are shown at the bot- 
tom of Fig. 01 Only a linear background has been sub- 
tracted. Comparing our single crystal data with the poly- 
crystalline results by Gordon et al. confirms the trivalent 
character of the cerium ion. We also find, similarly to 
Gordon et al., that the experimental dipole transitions 
are much broader than the higher multipole transitions. 
In Ref. 0, [l6| it was reasoned that this is due to inter- 
actions of the dipole final states with the continuum. 

The two measured directions in our single crystal ex- 
periment can be clearly distinguished. The experimental 
spectra agree very well with our theoretical predictions. 
In particular, the vector q dependence is in good agree- 
ment with the experiment. For example the inversion of 
the anisotropy at around 105 and 113 eV, as predicted in 
the calculations, is well reproduced. Here we note that for 
the calculation the Hartree-Fock values of the Slater inte- 
grals were reduced to about 68% for the 4/-4/ Coulomb 
interactions and to about 88% for the Ad-Af interactions 
to reproduce best the energy positions of the main spec- 
tral features. These values compare well to a 60% to 80% 
scaling of the atomic values needed to describe isotropic 
rare earth M^s (3d — > 4/) x-ray absorption spectra and 
account for configuration interaction effects not included 
in the Hartree-Fock scheme [3l], [32| ■ The simulated spec- 
tra were convoluted with a Lorentzian and a Gaussian 
with a FWHM=0.3 eV and 1.32 eV, respectively to ac- 
count for the intrinsic line width and experimental reso- 
lution. This experiment shows clearly that in CeCu2Si2 



a is negative so that the wings of the Ce 3+ ground-state 
orbital point in [110] direction. 

In conclusion, our simulation show that vector q de- 
pendent NIXS at the cerium N^s-edge in the limit of 
high |q| is sensitive to the crystal- field ground-state sym- 
metry. In addition, we show that the higher multipole 
character of the scattering at large momentum trans- 
fers gives insight into anisotropies which are higher than 
twofold so that the two possible orientations of a Tj or- 
bital in tetragonal symmetry can be distinguished. We 
have verified this experimentally by performing a NIXS 
experiment on CeCu2Si2 single crystals. The spectral line 
shape and the observed differences between the two mea- 
sured directions q||[100] and q||[110] are well explained 
by atomic full multiplet calculations. These unambigu- 
ously show that the ground state orbital in CeCu2Si2 is 
of the type T 7 = -|a|| ± 5/2) + VT~c7\ T 3/2). Thus, 
q dependent NIXS provides information which is not ac- 
cessible with dc susceptibility, neutron scattering and/or 
soft x-ray absorption spectroscopy and thereby completes 
the spectroscopic toolbox for determining the crystal- 
field ground-states in tetragonal cerium compounds. It 
is now possible and will be of great interest to investi- 
gate systematically the orientation of the ground-state 
orbitals of other tetragonal, heavy-fermion cerium com- 
pounds. Also, NIXS with hard x-rays has the advantage 
of a large penetration depth so that ultra high vacuum 
is obsolete, which allows the use of advanced sample en- 
vironments for pressure and/or low temperatures. 
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